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Coherent Structures 

By CHARLES G. SPEZIALE 1 

In order to develop more quantitative measures of coherent structures that would have compar- 
ative value over a range of experiments, it is essential that such measures be independent of the 
observer . It is only through such a general framework that theories with a fundamental predictive 
value can be developed. The triple decomposition 

<t> — ^ + <t>C + <t>T (1) 

(where is the mean, <f> c is the coherent part, and <j> r is the random part of any turbulent field <f>) 
serves this purpose. For a statistically steady turbulence which possesses coherent structures with 
a dominant temporal frequency / we can take (see Hussain 1983) 
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where <£' = <j> — f ; — if f , and < • > denotes the phase average. For more general turbulent flows 
(that are not necessarily statistically steady or do not possess coherent structures with a dominant 
temporal frequency), on can take 
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where an ensemble average is taken over N repeated experiments with the same initial and bound- 
ary conditions. The coherent part of the turbulence can be taken to be 


<j>c =< <t>*\E > 


(5) 


where < *| E > denotes a suitable conditional average of <j> 1 (i.e., an ensemble average over flow 
structures subject to the occurrence of some event E. With such triple decompositions, the 
coherent and random parts of the turbulence will be the same for all observers (see Speziale 1986). 

It should be noted that double decompositions (see Hussain 1983, 1986) give rise to coherent 
and random parts of any turbulent field <j> that, in general, depend on the observer. Double 
decompositions should therefore only be used when the mean flow vanishes or is negligibly small 
compared to the coherent notion (see Speziale 1986). Otherwise, one runs the risk of extracting 
flow structures that are overly biased by the observer. 

The equations of motion for the mean and coherent flow fields, based on the triple decomposition 
(1), can be written in the form (see Hussain 1983): 


Du{ 

~Dt 




( 6 ) 


I 


1 ICASE, NASA Langley Research Center 


i 



330 


C. G. Speziale 


DUc r, 

Dt 


dp c dui du Cl 

£7 + ,,v - *-i 5^ - *•> 


^d7 


dp 2 

h t/V 2 Hi 

OXi 


d 

dxj 


< Wfj Urj 


> 


( 7 ) 


where v is the kinematic viscosity of the fluid, p is the modified pressure, and u i. = u + u c + w r is 
the velocity field which is subject to the continuity equation which yields the constraints 


V u = 0 


( 8 ) 


V . u c = 0 (9) 

V • u r = 0 (10) 


In order to achieve closure of the equations of motion (6)- (10), the Reynolds stress terms 


< Ur t U rj >,U ri U Vj (11) 

need to be modeled. The time-averaged Reynolds stress u Ti u Tj can be modeled using the currently 
popular two-equation models or second-order closure models (see Launder, Reece and Rodi 1975 
and Lumley 1978). The phase-averaged Reynolds stress < u ri u T] > primarily serves as an energy 
drain on the coherent motion and thus it is plausible that it could be modeled using a gradient 
transport hypothesis (see Hussain 1983). Hence, eddy viscosity models of the form 


< u rt u Tj > — ~vt 
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can be considered where v? is an appropriate eddy viscosity (sufficiently far from solid bound- 
aries, the Smagorinsky model can be tried). This approach, which bears a certain resemblance 
to large-eddy simulations, has an advantage in that the coherent motion u c is calculated directly. 
Furthermore, the level of computation required is substantially less than that needed for a di- 
rect numerical simulation since a coarse mesh can be used (the fine-scale turbulence is modeled) 
and for some problems (e.g., turbulent mixing layers) the coherent motion is approximately two- 
dimensional. In my opinion, there is a good chance that such an approach could yield useful 
new information concerning the nature of coherent structures and it is well worth pursuing future 
investigations along these lines. 
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